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SUMMARY 



This study was carried out in an effort to determine the 
longitudinal stability derivatives of a Ryan Navion airplane using dynamic 
flight testing methods. This wa3 done by recording the transient response 
of the airplane to elevator forcing functions in the form of steps or 
pulses. Data reduction was carried out using the equations of motion 
method and the derivative method. Due to partial failure of the instru- 
mentation the results from the equations of motion method were disappointing. 
More successful results were obtained using the derivative method in con- 
junction with the useful data acquired from the equation of motion method. 

Concurrent with this report an investigation was carried out to 
find the stability derivatives of the same airplane from both theoretical 
considerations and from steady state flight testing. This investigation 
is described in Reference 7. An effort has been made to correlate the 
flight test results obtained from the dynamic and steady state methods and 
predictions made from theory. This correlation appears in Appendix B. 

Reasonable values of C L( ^ were obtained from dynamic testing 
which were consistent with speed, and position of center of gravity changes. 

Dynamic testing also gave reasonable values for the damping 
derivatives which were consistent at the sane speed, but which indicated 
an increase in the damping derivatives with increase in speed. 

The proper changes in values for the stick fixed maneuver margin 
were found as the position of the center of gravity was changed, however 
a decrease in maneuver margin was found as speed was increased. 

The elevator power was found with good consistency as the center 
of gravity was changed, but decreased in value as the speed was increased. 
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All values of Cm ^ were somewhat high as compared with theoretical values. 

An attempt was made to explain the changes in Cmj- and stick 
fixed maneuver margin that occurred with speed changes in terns of aircraft 
non-linearities. 
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LIST OF SYMBOLS 



G.* 

CmJ 



d£n 

2 



steady state lift coefficient. 

partial derivative of the moment coefficient with respect 
to angle of attack* 



n _ 

2d <. 



partial derivative of the moment coefficient with respect 
to = 'T' oc . 



C?r, -h-Qi partial derivative of the moment coefficient uith resoect 
d& W to de = 




c.- 



dC 

) K 



elevator power = partial derivative of the racnent 
coefficient with respect to elevator angle. 

slope of the lift curve 

(Pitching moment, ft. lbs.) 

'/z f d* S c " 



oC = perturbation angle of attack, radians. 

Q - perturbation angle of pitch (positive nose up) radians. 

do = ^ radians, 

d 0 - 'Z' & radians, 

7l - perturbation normal acceleration v positive for push over 
maneuvers ) . 

/* - §n d c /r = r^-fnc/d 
£ - perturbation elevator angle, radians. 



f~d = stick force, lbs. 

rate of change of downwash with angle of attack. 
V - velocity, feet per second. 



v 



s/c 


mass of the airplane, slugs. 


OC - 


«/ _ 
//5C 



t' - seconds. 





2 

moment of inertia of the airplane about the Y-axis, slug ft. 


A - 


ZK?/ z* 



ky r y -A At , radius of gyration, feet. 



c -- 


mean aerodynamic chord, feet. 


5 - 


wing area, square feet. 


i- 


tail length, distance from the center of gravity to the aero- 
dynamic center of the tail, feet. 


Sw~ 


2 

134.2 ft. = wing area. 


St - 


2 

43 ft = horizontal tail area. 


Se- 


2 

15.04 ft * elevator area. 


A w 


6.04 = aspect ratio of wing 


/U- 


3.98 = aspect ratio of horizontal tail. 


H 

> 


33.33 ft = wing span 


4^ = 


13.17 ft = tail span 


C -- 

O IV 


5.7 ft - mean aerodynamic chord of wing. 




3.34 ft = mean aerodynamic chord of horizontal tail. 


f - 


15.04 ft = tail length, distance between center of gravity and 
aerodynamic center of the tail. 


^ 0- 


■ .24 - position of wing aerodynamic center in % 


1 iv - 


2° at root = incidence of wing 
1° at tip 


4 - 


-3° = incidence of horizontal stabilizer. 



Dihedral = 7.5° 
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INTRODUCTION 



Dynamic testing of aircraft for stability parameters has received 
increased emphasis with the advent of guided missiles, sonic, and super- 
sonic aircraft. Some of the major efforts in this field have been dynamic 
testing of a 3-25 by Marshall E. Mullins at Edwards Air Force Base, (Ref. 

4)i the dynamic testing of an F-30A by R. C. Kidder at Cornell Aeronautical 
Laboratories, Inc. (Ref. 5)* and guided missile testing by The Naval Air 
Missile Test Center, Pt. Mugu, California. In all previous cases no effort 
has been made to compare the results with those obtained by steady state 
flight test methods for the same aircraft. This report is part of a 
triple project on a Ryan Navion in which dynamic results are to be compared 
with theoretical and steady state results. The actual comparisons are 
included in Apoendix B. 

The dynamic testing with which this report is concerned consists 
of analysis of the transient response of the aircraft to an elevator forcing 
function. Using various forms of the equations of motion the stability 
parameters for the airplane were reduced from the data by the method of 
averages, least squares, and by the use of graphical plots. 
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THEOHY AND ANALYSIS 

The normal procedure in the determination of aircraft transient 
motions is to solve the differential equations which describe the aircraft 
motion. In doing this the constant coefficients in the equations of motion 
are estimated from theoretical considerations or measured by steady state 
test methods. 



stated above: ie., the determination of the coefficients in the equations 
of motion if the time history of the aircraft transient motions are known. 
The response history of the aircraft can oe displayed as a function of 
frequency as well as time. In theory the response in either the time or 
frequency domain can be analyzed to determine the constant coefficients 
in the equations of motion. This report deals with the determination of 
the constant coefficients of the longitudinal equations of motion from 
transient responses displayed in the time domain. 

Development of Equations: 

Throughout this report it is assumed that the incremental changes 
in forward velocity are negligible. With this assumption the equations 
of motion as developed in Ref. 1 are: 



The problem dealt with in this repcrt is the inverse of that 



o{ + d - dO - 0 



( 1 ) 




(2) 
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If elevator lift is considered important, its effect can be included in 
the lift equation since: 

- Cmj 






4, 

'C 



( 5 ) 



Substitution in equation (1) yields: 



cx'+dc< - de - £zL 

Z 






o 



(4) 



Solving equation (4) for cJo( and substitution in equation (2) yields: 




Normal acceleration, n, and elevator angle, S , may be substituted for <X 
since: 



<*- - S-n + Cm 



C 



L oC 






l s 



Substitution of this expression for 0< into equation (5) yields: 



( 6 ) 




C 



"'du. 



n + ( C "‘«l» + * C ”i + - 0 



(7) 
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This equation if integrated once with respect to 



yields: 




It is possible also to eliminate all variables except n, 3 and their 
derivatives from equations (2) and (4). Solving equation (4) for 39 and 
(3 l 8 and substituting these into equation (2) yields: 



(c~„+ 









' < - m s (■ ' ' 
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_ hCiu 


)c/o(' - Ac/ ex' 




■‘)s 


h Cm £ 
+ 3^/c 


o! £ = 


£> 


(9) 


for 


d o( , and 


do( 


in terns of n and £ yields: 




o)o( = 


■ -A- cJn + 

C Lo ( 


C L ^/ c 


dj 


(10) 


d* 


- -Cl 
Cl, 


cTn - 


f. J 

Cl, 4 


dV 


(11) 



Substituting equations (6), (10), and (11) into equation (9) yields: 



dVl + Cm^ ( 






K A/ c C L h Cl 
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-%■ hC L l 
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Equation (12) if integrated once may be expressed as: 




where 



$ -/ h d H 

(14) 

Throughout the remainder of this report the following quantities 



IT - ( 



Ijn- JS + Cr (r 
' c L p c crkv"^ ' 



are defined as: 
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A = 


f' / C rT) ot ^ ol<x) 

' * J 


(15) 


B = 




(16) 


C = 


C "‘( " £*) 


(17) 


L = 


<4 ( + c ye) 

c B a/ n j 


(13) 


M = 


r' / C/_ o( , Cm oc ^ 
Cm j( C L Q 4 ) 


(19) 



Substituting equation (15) through (19) into equations (7), (3) and (13) 
yields: 

An -t-ScIO + CS-hcJ^ ( 20 ) 

Afn + Be + Cfs - Bc/e (21) 

Lf n +/y /p =y U2) 



Equation (22) may be further rearranged as: 



L 



ir 



+ 









IT 



= / 



(23) 



Solution of Equations: 

In order to solve the equations of motion for values of the deriva- 
tives, one quantity must be known. In this case the known quantity was the 
airplanes moment of inertia about the Y axis. To determine the moment of 
inertia, the airplane was oscillated about the jack points as shown on Fig. 1. 
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If it is assumed that the airplane as suspended represents a second order 
system with zero damping, it is possible to determine the moment of inertia 
about the jack points, from the formula shown on Fig. 1. The moment of 
inertia was then transferred to the center of gravity. The position of 
the center of gravity was determined by using scales under each jack point 
and under the tail for two different inclinations of the airplane. With 

the moment of inertia about the center of gravity known it was possible to 

/ 

determine 4 about the center of gravity. By keeping an accurate record 

of weight added to and removed from the airplane A was determined for each 

run. The accurate determination of this quantity for each run is very 

important as small errors in A cause large errors in the solution for the 

constants in the equations of motion, equations (20), (21), or (23). 

Examination of the equations of motion 3hows it would be most 

desirable to measure £e , , dd , and d*Q directly. Since an angular 

2 

accelerometer was not available, direct measurement of d 0 was impossible. 

Equation (6) shows the desirability for having normal acceleration data in 

order to determine Ci^ . For this reason, and realizing that there are 

certain inherent difficulties associated with the accurate determination of 

a, it was decided to use normal acceleration as one of the measured 

variables. In order to make the same instrumentation system available for 

use with steady state data and so that the elevator hinge moment equation 

might be considered in future studies, it was decided to include stick 

force as one of the measured variables. The resulting instrumentation 

yielded simultaneous values of a, n, A e , d0, and F . The actual circuits 

s 

used are described in Appendix A. 
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Since h is a measured quantity, measurement of n, , d9, and 

2 

d 9 at three different values of time would permit the solution of equation 



at three values of time would permit the solution of equation (21) for 



would permit the solution of equation (22) for L and M. 

Methods of handling redundant datai 

If the data were perfect, evaluation of the variables at only 
three instants of time would be sufficient to uniquely determine the 
coefficients of equations (20) and (21). With less than perfect data, 
however, some method must be used to determine the set of coefficients which 
most logically fit the data. The authors have used several methods of 
handling the redundancy arising from evaluation of the variables at many 
values of time. The first was a simple application of least squares assuming 
the error to be all in one variable i the second was a method of averages i 
the third was a graphical plotting method. This latter method is the most 
desirable since it gives a physical picture of the data scatter and the 
fairing process. The lack of physical feeling for the fairing process 
when least squares or averages are used to reduce the data is of great 
importance . 

When the method of least squares was used the normal equations of 

/ \ 2 
equation (20) assuming all errors in d 9 ares 



(20). for A, B, and C. In like manner measurement of J n, 9, and d9 



A, B, and C. Knowledge of 




j , and v for two values of time 



A w + B £}'/( dd + C B,7l£ - A d A) 

A ZncJ9 + +CZ&d&- A 

A S +- B<CSd9 ¥- C i ~ dl9 



(24) 



(25) 

(26) 
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In applying the method of least squares it was necessary to carry as many 
places as possible in combining the variables and solving the simultaneous 
equations. No justification is offered for the apparent inconsistency 
of carrying eight to ten significant figures in the summations of the 
various combinations of variables when the variables themselves cannot be 
taken from data sheets with more than three significant figures. It is 
merely stated as a finding that as many places as possible must be used 
in order to get reasonable answers using least squares procedures. It 
follows that least squares procedures should not be attempted without the 
use of an eight t6 ten place calculating machine. It is also to be noted 
as a corollary of the above that least squares procedures are very unfor- 
giving of small errors in data transcription, summations, etc. 

The method of averages consists of simply using the average values 
2 

of n, d9, £ , and d 0 taken from random groupings of the Individual data. 

The best results were obtained by simply grouping the data into approximately 
even thirds. The first of the three simultaneous equations then would be 
made up of the summations of the variables over the first third of the time 
traces considered appropriate. This method consisted of dividing the data 
into three such groups. The simplicity and speed are such great advantages 
over even the simplest application of least squares that its applicability 
should be at least evaluated in any extensive data reduction problems. 

Practical Considerations 

Data for eight runs were reduced using any or several of the 
methods described above. The results were consistent for different methods 
applied to the same run but the only answer which seemed appropriate was 
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the value for B = + Cn^ a ). The values of A and C were of improper 

magnitude and inconsistent from run to run. With this in mind the equations 
for the unknowns A, B, C in determinant form were examined. 
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Examination of equations (27), (23), and (29) shows that equation 
(23) is the only one where an error in the pitch rate gyro calibration 
would not affect the results. 

An examination of the pitch rate gyro revealed that considerable 
friction had developed in the bearings and it was very doubtful if the 
gyro wheel was getting up to governing speed. 

If it is assumed that the gyro maintains a constant speed during 
the 1.5 seconds during which data is being taken for any particular run, 
the reduced values of B will be correct. Values of A and C will of course 

be directly dependent on the gyro speed for that particular run. 

Because of time limitations it was impossible to repair the gyro, 
recalibrate, and take new data. This necessitated re-exanination of the 
equations of motion in order to eliminate pitch rate as an unknown. The 
method is very similar to that proposed in Ref. 2 and called the derivative 
method . 

An examination of equation (14) shows that if (C //?cyx + ) and 

C* are known, an estimation of C 77Ii will define the quantity if • 

Equation (23) is then recognized as that of a straight line where the 
variables are ^ n /if and ^/if and the intercepts with the coordinate axes 
determine values of L and M. At each instant of time it is only necessary 

to compute f? c , f£ , and if \ make a plot of equation (13) on linear 
coordinates! and fair the best straight line through the points to determine 
the intercepts. These graphical plots were extremely enlightening. They 
indicated that L can be determined to much greater accuracy than can M. 

A blind application of least squares or averages could never indicate this 
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condition. Knowledge of this condition would be valuable for deciding / n 
which variable to minimize the error for application of least squares. 

It should also be noted that the answers are very insensitive to 
the estimated values used since the known value of (C^ + makes up 

the major part of v. 

In the process of data reduction it was attempted, with poor 
results, to use a data instead of n for solving the moment equation. The 
fifteen hundred series runs yielded very poor results. This was accounted 
for by the fact that an accurate in-flight calibration of the angle of 
attack indicator to account for vane position error effects was never 
attained for these runs. The nature of the difficulty is described in 
Appendix A. 

The a data was used in equation (6) to obtain values of . The 

a 

use of a data necessitates taking into account the effect of up wash and 
the effect of pitching velocity. The effect of up wash is determined from 
steady flight calibration of the angle of attack indicator. The effect 
of pitching velocity is determined by measuring the distance of the vane 
ahead of the center of gravity and computing the velocity component at the 
vane due to pitch rate. Another effect that it was not possible to account 
for quantitatively was the effect of acceleration in pitch. Due to flex- 
ibility of the boon angular acceleration of the airplane caused bending of 

the boom in such a manner that individual values of the C. determined at 

a 

each instant of time yielded a plot as shown below. Better stiffening of 
the boom should eliminate this source of error. 

As the dynamic response of the angle of attack vane was not deter- 
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mined, any natural dynamic response of the vane consisted of an error in 

the a data of undetermined magnitude. 

For the above reasons a data was used only to determine from 

a 

the sixteen hundred runs. 







The solution of the equations of motion by any of the techniques 
used does not allow solution for unique values of the derivatives Cyn do c , 
Cm^ » * and • In order that the results may contain 

approximate values for the above derivatives it was assumed that 
and = .5. It is important to realize that without the above assumption 

Ofot 

dynamic flight test procedures as described herein do not yield unique 
answers for the individual derivatives. 
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PROCEDURE 

The procedure followed in this report may be considered in the 
following parts: 1. Instrumentation, 2. Calibration, 3. Flight test, 

4. Data reduction. 

1* Instrumentation 

The design, construction, and installation of the instruments 
used in this report took approximately four months. The circuits actually 
used as well as the design criteria may be found in Appendix I. The 
reasons for the choice of variables to be measured may be found in the 
Theory and Analysis section of this report. 

2. Calibration 

The calibration curves for each circuit used may also be found in 
Appendix A. A brief explanation of inherent difficulties in measuring 
certain variables may be found in the Theory and Analysis section of this 
report. Calibration was, of course, carried out continuously as the flight 
test progressed to insure that the calibration curves remained applicable. 
The one exception to the above statement was the calibration of the pitch 
rate gyro. It was calibrated in the final stages of instrument installation 
at the Naval Air Test Center, Patuxent River, Maryland, and no further 
calibrations were performed. The difficulties encountered in the early 
stages of data reduction emphasized the desirability of continuous calibra- 
tion checking. 

3. Flight test 

The flight procedure for transient dynamic testing is extremely 
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simple and represents the greatest advantage of dynamic testing over 
steady state testing methods. Smooth air is essential in order that a 
good steady state flight condition may be attained. After attainment of 
a good steady state condition, the steady state is disturbed in some pre- 
determined manner. Pulse type elevator forcing functions were applied by 
simply moving the yoke a predetermined amount, holding it in the new posi- 
tion for approximately one second, and then returning it to the original 
position. The resulting pulses were poor approximations to an analytical 
pulse. However, in an equations of motion technique this was not important. 
The best method found for applying a step elevator input consisted of 
holding the stick against the force caused by a preset value of trim tab 
and then simply letting go of the stick. The elevator is so heavily 
damped that the resulting elevator motion was very nearly a step. 

All data was taken at 19.5 inches manifold pressure and about 1850 
engine revolutions per minute. Low engine speeds were required to reduce 
angle of attack boom resonance. Variations in speed were obtained by 
varying the rate of descent. Data was taken for two center of gravity 
positions and two speeds using both pulse and step elevator forcing functions. 
4. Data reduction 

The various methods of data reduction are described in the Theory 
and Analysis section of this report. 
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SAMPLE CALCULATIONS 

A sample of each type of calculation is carried out in this section 
to further clarify the procedure used. The calculations are carried out 
on run #1640. The values of the variables taken from run #1640 may be 
found in Table II. Table I contains the various constants required in the 
calculations and the calibration factors taken from the calibration curves. 
Equation (21) is 

V ^ + £ e + cf< J = A cJ& 

The values for Jn « 0 , JS anJJ6 may be obtained for each point indicated 
in Table II since; 

Jn di /r = rfvcJt - Jr J nc/ t scoonc/s) 

where C ^ ~ calibration factor for n in 

Jd c/ /p - — L JSc/f ( inch seconds) 

fdoS = # = C 6e fee// A, \r>ch Seco nds) 

de = 5 % = - T C ds 9 fine hts) 



With values of the variables at each point three simultaneous 
equations may be obtained as follows where the summations are over the 
various points of time; 
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A'zfn’*^. * B Z e 
A±fndi- r * 3 ie f 

S S 

Alfnc/^zr + £ ze + 

9 9 



Cd J- J cf d /jr = Jl £ do 
C i /JV*r - hkde 
ckj^j^z - A £ do 



Actual calculations yielded the following simultaneous equations; 



- .000091 A + .017371 B - .020173 C = .010376 

- .198996 A + .110355 B - .074817 C = .020030 

- .601369 A + .204940 B - .125027 C = .014376 



Solution of these equations simultaneously for B gives B = - .141341. 

It is clear that equation (7) as it stands could be solved using 
2 

values of £ , n, d©, and d ©. These variables in inches and seconds units 
are also included in Table II. 

If equation (6) is rearranged it is possible to solve for . 

■ r . - C n + € 

" g 7c 

Of 



Values of n, S , and o< may be obtained from Table II using the calibration 
factors of Table I. The value of 0C so obtained should be further cor- 
rected for the error introduced in the vane indicator by pitching velocity. 




/?<?/£ //U<? A V/nrJ 



y* s ft 



A of 
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where V is true velocity in ft/sec. 

It is necessary' to estimate C& ^ was estimated theoretically to 

be -1.3. An incorrect value of C //t g will of course be reflected in an 

incorrect value of but the percentage error in will be much less 

a a 

than the percentage error in C „ . ^ . Using all points of Table II the 

following equation was obtained: 

n ~ C c 71 f Cms '<? c 

' ^/C r h 

‘4* 

^ PC 

i 

C. - -( -7+ ?) (-3.4-4-79)+ 

L «' 7 Z.Q4- 

. S3 7 




7?. (c> C> -f . /£, 
■ 7T37 



s.oc 



For this particular case an increase in of 50$ gives an 

increase of 12$ in . 

“a 

The terras on the right hand side of equation (13) have been esti- 
mated as shown at the end of this section for runs #1522, 1523, 1640, and 
1652. These estimated quantities have been determined using values of 
( Cw + ) determined by equations of motion method and as rationalized 

d <?C 

in Fig. 2. An average value of C/ ^ and the theoretical value of Ow$~~/3 
have also been used in these estimations. Tables III, IV, V, and VI contain 
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the calculations for v and the variables 



and 




These variables 



are plotted and the values of L and M determined on Figs. 3 and 4. 
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RESULTS 

The values of C. were consistent with theoretical considerations 
<* 

and no trends were indicated as speed and center of gravity position were 
changed. 

The damping derivatives were found with good consistency at the 

same speed but an increase in damping derivatives ( C ni + Cfy, ) was found 

ax dy 

with increase in speed. 

r C ni Cl 

The values of stick fixed maneuver margin £ — -?• + 

indicated a change of the proper amount when center of gravity position 

was changed at constant speed. However, a decrease in stick fixed maneuver 

margin was found as speed was increased. The stick fixed maneuver point 

for 90 mph was 58$, for 110 mph it was 52$. 

The value of was with some consistency at each speed, 

& _ 

but a decrease was indicated in as speed was increased. The values 

of C„ were somewhat high in comparison with theoretical values. 

i 

The only immediately apparent enigmas in the results then seemed 

to arise as a result of speed changes. This would lead one to suspect an 

error in some or all of the constants used in data reduction which were 

functions of speed. These constants are C * an<^ ^ ie va l ues °f 2~ an( ^ 

C used were very carefully examined and no errors could be found. 

Yj 

During the in-flight calibration of the angle of attack indicator 
a condition was found which might explain some of these changes with speed. 

A sharp increase in up wash at the wing tip was noted as speed was decreased 
below 85 mph indicated. The test aircraft was equipped with a spoiler near 
the wing root. The possibility existed, then, that non-linearities in flow 






I 

I 
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and lift were occurring during pull ups froi the 90 mph steady state. 

If non-linearities in the flow did exist, they would be reflected 
in low values of normal acceleration. Further the actual error in the 
normal acceleration would become greater with time. The recorded values 
of normal acceleration would be less by a progressively larger amount than 
that which would have been recorded sans non-linearities. The recorded 
values of would have been less but not by as large a percentage as 

normal acceleration itself. The error in was not as apparent as that 

for n and J ' n . 

The equation from which and stick fixed maneuver margin were 

obtained was: 



/ 



Sl 






S 

Sa 



= / 



, /■ M ^ 

ly i/ 

The quantity v as determined in Tables III, IV, V, and VI was predominantly 
made up of n, particularly after the first few tenths of a second. The 






V 

was a 



non-linearities discussed then, if present, would have tended to make 

and & too large, the error increasing with time. Since JA 
V Cr -V 

larger quantity than these errors would tend to decrease the slope 

'V 

of the lines faired into the data in Figs. 3 and 4. This could have given 
the increased apparent stability found at lower speeds. Also since the 
points at later times were used predominantly in fairing Figs. 3 and 4, 
the increased value of at 90 mph could also be caused by these non- 



linearities. 



An attempt to evaluate the effects of non-linearities on the values 
of the damping derivatives as determined by equations of motion method lead 
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to no definite conclusions. It is possible that 3one of the apparent 
increase in damping derivatives with speed may have been caused by low 
values of normal acceleration at the lower speed. The increase in damping 
derivatives with speed indicated an increase in tail efficiency. The test 
procedure required as low an RPM as possible to reduce the angle of attack 
boom resonance. Power used was 19.5 inches, full low RR4 at both speeds. 
This required an increase in rate of descent with increase in speed. This 
does not appear consistent with the increase in tail efficiency unless the 
tail wa3 operating in a particularly unfortunate position with respect to 
wing wake at the lower speed. 
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CONCLUSIONS AND RECCKMENDAT IONS 

Itynamic flight testing is in general a perfect ly feasible method 
of accurately determining the stability derivatives of an airplane. 

Certain basic requirements for obtaining good results are as follows: 

1. Inclusion of pitch rate as one of the measured variables. 

2. Weight, moment of inertia, and airspeed must be known very 
accurately. 

3. Instrumentation must be very accurate without unknown instru- 
ment dynamic characteristics. 

Dynamic testing has the following advantages and disadvantages. 
Advantages: 

1. Simplicity of flight procedure. 

2. Short duration of in flight testing. 

3. Non-linear effects may be included in the equations of motion. 

4. Values of the derivatives obtained are the ones the airplane 
feels in response to transient maneuvers. They are the values for velocity 
equal a constant. 

5. Close correlation exists between frequency response determined 
from dynamic testing and auto pilot design requirements. 

Disadvantages: 

1. Instrumentation is generally very expensive. 

2. Data reduction is long, tedious, and subject to error. 

3. Some of the results can not be compared to static flight test 

2 

results which are for C^V equal a constant. 

In recommending that the work be continued the most important 
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improvements in the methods used are listed below, 

1. Further efforts should be made to obtain accurate pitch rate 

data. 



2. A pendulum or turntable should be constructed for calibration 
of the rate gyro and accelerometer. 

2 

3. An angular accelerometer would permit the measurement of d 0 

directly. 

4. Angle of attach data would be improved by boom stiffening, 
and dynamic calibration of the vane. 

5. A frequency analysis of the airplane response might yield 
interesting results. 
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APPSNDK A 

Description of Test Equipment 

The test airplane was a standard Ryan Navion instrumented to record 
time histories of elevator angle, angle of attack, stick force, pitch rate 
and normal acceleration. These variables were recorded with a Consolidated 
Engineering Corporation type 5-116-P3-14 recording oscillograph. A wiring 
diagram for the instrumentation system Vised is shown in Fig. 5. Schematic 
diagrams showing the operation of the individual circuits appear as Figs. 

6 and 7 . 

Certain basic considerations were common to the design and operation 
of all of the circuits. These were: 

(a) To make the equivalent resistance of each circuit be 350 ohms 
as seen from the oscillograph leads, in order to meet the manufacturer's 
recommendation for 62% of critical damping for optimum galvanometer response. 

(b) To operate the sensitivity adjustment of each circuit at 
about 50fe attenuation. 

(c) To use as small balancing resistors as possible in order to 
maintain linearity of the output, and at the same time keep the resistors 
large enough to limit the current drawn from the battery. 

(d) To obtain maximum galvanometer deflection for the range of 
variables measured at each flight condition. The galvanometer deflection 
was one inch per 12.4 fJ. amps of current. 

Power for the instruments was obtained from a 45 volt dry cell battery 
with a center tap. Voltage from this battery to the bus lines was regulated 
to 20 volts for all tests using a variable resistor and voltmeter. 



-30- 



Slevator A ngle and Angl e of A ttack Circuits . The ^ and o C circuits 
are shown schematically in Fig. 7. The potentiometers used to sense the 
angular motions were 12,500 ohm, Micro-Max precision potentiometers with 
outputs guaranteed to be linear within one per cent. The motion of the <x 
vane and the elevator was transmitted to the potentiometer shafts in a one 
to one ratio. The 2500 ohm balancing potentiometer is used to adjust the 
circuit output to zero at various elevator or a-vane positions. The addi- 
tional 5000 ohm resistors in series with the balancing potentiometer were 
required to limit the current drawn from the dry cell. The 1500 ohm load 
resistor was required to adjust the current to the galvanometer to the 
desired magnitude. Sensitivity adjustment was achieved through a 500 ohm 
impedance pad. The 1000 ohm variable shunt resistor was placed in parallel 
with the circuit to make the galvanometer external damping resistance be 
350 ohms as required. 

The calibration of the elevator angle circuit is shown- in Figs. 8 
and 9. These calibrations were made using a propeller protractor mounted 
on the elevator to measure A d c s . Fine adjustments of the elevator 
position were made by moving the elevator with a hydraulic jack under the 
trailing edge. 

Calibration of the angle of attack circuit was made in flight using 
the propeller protractor to measure A . Calibrations were also made 
on the ground using a specially made fitting on the boom to measure the 
angle between the vane position and the boom centerline. These calibrations 
are shown in Figs. 10 and 11, and illustrate the difference in calibration 
under flight conditions caused by airloads on the vane and boom. Due to 
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an oscillograph breakdown it was impossible to get an in-flight calibration 
for the 1500 series runs before the sensitivity and zero balance settings 
were changed. Therefore the in-flight calibration for the 1500 series runs 
was made using the difference between ground and flight calibrations for 
the 1600 series runs and extrapolating for a correction to be added to the 
1500 series ground calibration. 

S tick Force . Pitch R ate and Normal A cceleratio n C ircuits . The oper- 
ation of the F g , 9. and n circuits are similar and are shown schematically 
in Fig. 6. The strain gauge bridge circuits used to measure these variables 
will be described individually later. The 5000 or 50,000 ohm balancing 
potentiometer in each circuit gave a zero output adjustment to the circuit 

m 

over a wide range of values of F ,0, and n. The 5000 or 50,000 ohm re- 

3 

sistors in series with the balancing potentiometers were necessary to provide 
a current limit to protect the galvanometer in case the circuit should 
inadvertently be turned on with the sensitivity adjustment full open and 
the balancing potentiometer turned against the stops in either direction. 

The sensitivity adjustment and damping resistors are the same as described 
for the a and circuits. 

The accelerometer used in the n circuit was a model F-l. 5-380 accel- 
erometer manufactured by the Statham Laboratories, Incorporated. The manu- 
facturer' s specifications required an 13 volt input to the accelerometer 
which was obtained with a 100 ohm variable resistor in series with the 
circuit leads to the power supply. The accelerometer had a range of +1.5 g' s. 
The circuit was calibrated between zero and one g and assumed to be linear 
to 1.5 g. The calibration curve is shown in Fig. 12. The noise level in 
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the n circuit was extremely sensitive to airplane vibrations, particularly 
to engine rpra. In order to minimize this the 7r-section filter shown in 
Fig. 6 was used to filter out vibrations above 12 cycles per second. 

The stick force bridge was made up of 3R-4 strain gauges cemented 
to the spokes of a specially constructed force wheel which replaced the 
standard pilots wheel in the airplane. Calibration of the stick force 
circuit is shown in Fig. 13. 

The pitch rate circuit bridge was made up of SR-4 strain gauges 
mounted on a leaf spring which restrained the precession of a gyroscope 
when the gyroscope was subjected to an angular velocity. The SR-4 strain 
gauge bridge requires a 15 volt input which was obtained by using a 100 ohm 
variable resistor in series with the circuit leads to the power supply. 

The pitch rate gyro was calibrated on a calibrating pendulum at the Naval 
Air Test Center, Patuxent River, Md. , and subsequent recalibrations were 
not feasible. This calibration established the linearity of the circuit, 
and is shown in Fig. 14. This linearity extended through six points not 
appearing within the range of Fig. 14. A 250,000 ohm +0.1$ and a 50,000 ohm 
+0.1$ precision resistor was placed in parallel with one side of the bridge 
as shown in Figs. 5 and 6. At the tine of the original pendulum calibra- 
tion the bridge unbalance caused by these calibration resistors was noted 
to give a circuit output equivalent to particular values of pitch rate. 

Then with the instrument installed in the airplane a calibration could be 
made for each test flight. This was done by momentarily closing a calibra- 
tion resistor switch with the aircraft in steady flight. A galvanometer 
deflection was thus obtained for a known value of pitch rate. This point 
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along with the point of zero deflection for zero pitch rate was then used 
to establish a calibration curve for the flight. 

The electric motor in the gyroscope was powered from a 24 volt 
storage battery. It was discovered during the data reduction that the 
gyro motor had had a progressive bearing failure, and that 0 values were 
unreliable. This failure could have been discovered had a calibrating 
pendulum been available for calibration after each test flight. 
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APPENDIX A 

The following is a tabulation of the results from 
Princeton Reports No. 231 and 232. Report No. 232 consists 
of the determination of longitudinal stability parameters, 
for a standard Ryan Navion, using theoretical calculations 
and steady state flight test- techniques. Report No. 231 
consists of determining as many as possible of these same 
parameter s for the same airplane using dynamic flight test- 
ing techniques. 
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-0.031 


N o 


0.380 
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0.440 
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N m 
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0.52 
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Clean - Power Off 
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.xonlnatlon of the t bloc *'bovo ohovrr tho folio**iny: 

1 . Taluos of C/. oter ined bv the different methods 
ore 'ill -rithln 1 °'. 

£• Ti*eor .tic 0. °nd to vj t t A v^lucc of k^ngra© 

vsry closely vhtle tho v~l ”0 of fro*. dyn” 0 ’ic to;tiny coe c 

to bo ?''ont 30 "' hijh. doo“ of this could be " 0 counted for 

b” t u o fret thrt no '-lire'' ritiee in tho elevator enylo 

po to 'tiorotor used for '’no ic verb cause' •* sc.scibi 1 it** 

of 1 ''-I'"- ' error in iuividual v°lu®c of f’<^ 6 lovctor on **lc . 

/s 

'o*'o* r or, it^isi robable that the clo e of t -,o 0 -lo , ” , tor nr, .';lo 
c^libr-" tion curvo c^uld bo core thrn 10" off. In steady 
st- to testing, .-o.*er for level flight *”*c uccd. In the 
b*n n nlc testing, lo ,s co*.*cr v r o s used in order to ’coop tho 
c”> ;ino T lov to ’'.void -njlo of ’'tt’ch boon rocor-once. 

The 110 ’ d sir oeeod -.*-uc attained by est”bliohin 3 0 f *to'*dy 
r'-'.to of descent. 

3. T’-'O voluec of the ctlch fi od maneuver oint //*> 
ny.air. prec bet-noon theoretical and steady state values ’>*hile 
t' e v-alue f ran d* r nauic testing ceenc t r be 17 hiyh. The 
difforeroe rooms. e: cor-ivo since it ro^ro cents n 47 ' 'lescrep- 
■* , nc*r in the •’’cneuver r°ryin. The yrafaical lots of Fiy**. 

3 and 4 of ,:, o ort o. °31 indicate the difficulty of obt'-ln- 
iny fro - ’ 1 lyn-nic d°ta the intercept "hich determine t N e o~n- 
^uvor s’T :in. >en th *»u *b tho inter'ce t i defined, by lines 
cro^ciny at ~ c*a.l on pie this c~n not '"’count for t'~e 47* 
disore. •”ncy . Oy n an i c - or’: * r i c .veroly , m**’ -''r^d by failure 
of t l e r-’te yyro -md continued invo itiy tlon *dth cc^'^te 
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rate yyro measurements may yield better res 'It . The efloct 
of c-’rrylny tho lo .•or cover for the dyr^nic ton tiny is °1 .o 
apearont when tho povor-off maneuver aryin iroa steady state 
results is considered. If tho values from dynamic testiny 
are compared with the po-’er-off values of steady "t°te results 
tho di- srepc.ncy is °7 ' indicatiny that the effect of carrying 
tho lo' r er o-:er In dynamic testiny Day he very l°ryo. 

The vlue of oht-inod fron steady st°te testin ' 
is rue-'' lo ’or th°n t'^-t from dynamic testiny for the sane 
re'- sons cousin p tho discrepancy in tie maneuver m~ryin. 

Tho difference in Cm ^ is further n ccontu°ted hy the 

snail valu^ of t s o d^mpiny derivatives obtained from steady 
3 tate too tiny. Tho thoorotio-l value of Cm*. determined 
from the stic’.: fired neutral ,oint ay re os with the v^lue 
of Cn^ from ste-dy state te-stiny and the theoretical 
value of Cm^ determined from otic-: fir.od maneuver point 
ayroes fnily veil with the value of Cm * determined fr~~om 
dynamic tests. This show's that dynamic test results are 
very closely related to maneuverin' - ; stability while the 
stench” "t'' 4 ’e results are more nearly rol t^d to static 
stability. 

The values of the dam piny derivatives ana 

abtained from theory and dynamic tostiny ayree closely while 
those obtained from steady state tr-tiny •’re too lo - ’. This 
con ho primarily accounted for by the f°ct th-t tho d-np- 
iny derivatives are very diffic lit to e.-tr-ct from the 



steady st n to test data and ■'.re proswed to- be in error 
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4. T'-e TO"t a 't v~ri"nce "^Gt’ eer t’-o th.or ti” 1 and 
steady state flight toot rosults, ---ith •: ho e- ©option of 
parameters rovi • usly discussed, -’ns found to be in the 
stick free maneuver point for t '10 go -er-on, clean condi- 
tion, ~nd in th© stick flood neutral oint for tho jo’-or- 
off , l^ndin y condition . 

It is b„liev.3d that tho discrepancy in th© ..tick fro© 
maneuver joint any do attributed to th© f->ot tlvt at lo” 
nor "l "c colorations tho stick fore© required i° lo-r and 
erroneous ind _c?.tions of stick ^orce are .ossicle because 
of the friction in the control column and elevator system. 
Oo.tr. reduction and analysis indicated no change in naneu- 
vsriny oint betveen nor r ' , ' , l accelerations of 1.0 to 1.3. 
^ouovor, if an ' ntorool^tion is v- he between tho '-tick 
free maneuver joint for noml accelor- tions of 1.0 to 
1.5» " value cf -• nouvor point of ^ ^roaim^tely 33 m can 
no ro ay n noiic chord is obtained which corresponds ojaite cIogo 
with the 33 mean -".oro dynamic chord obtained Cron theoret- 
ical considerations. 

Tho stick fiacd neutral joint for the ouer-off, land- 
iny condition which was obtained from steady st^t<? to 'ts 
corresponds to t’~o trends indicated in t >e other jo- r er con- 
oitio^ f ' and configurations. 'kruever, the theoretical n r,s- 
lysis inlic-’tod conoid r^bly norc stability. ‘ s'-lynmont 
of t‘-o error involved can not ’x? don o with ~ny a^yroe of 

^ cc”-^' cy : ho- f ©vcr , it i3 believed tho nn'or ^r^i^n of + ho 

AAe 

diffe: cnc can bo ■ coo irtod for in the gt lu-ti^n of^f actor 
( / - - d Ajj . } la tho tho ore ti cal '•nod -ni't, 
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Unreliable . See description of test 
equipment section of report 
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